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ABSTRACT 


Inelastic  electron  scattering  experiments  were  carried 

7 
out  on  Li  within  a  range  of  momentum  transfer  squared  of 

2        -2 
.2<q  <  .65(F   ).   Inelastic  cross  sections  and  form  factors 

were  determined  for  the  4.55  MeV  inelastic  peak  and  the  in- 
elastic continuum  caused  by  the  disintegration  of  Li  into 
an  alpha  particle  and  a  triton  above  the  threshold  of  2.45 
MeV. 
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I.   INTRODUCTION 

6         7 
Inelastic  electron  scattering  experiments  on   Li  and   Li 

have  been  conducted  by  Bernheim  and  Bishop  [Ref.  1,  2]  in 
which  the  inelastic  form  factors  were  determined  for  these 
isotopes.   These  same  authors  [Ref.  3],  as  well  as  Hutcheon 
and  Caplan  [Ref.  4],  have  also  investigated  the  inelastic 
form  factors  of  the  continuum  at  excitation  energies  above 
the  a  -  triton  threshold  for   Li  and  above  the  a  -  deutron 
threshold  for   Li.   Bergstrom  and  Tomusiak  [Ref.  5]  of  the 
Saskatchewan  Accelerator  Laboratory  presented  a  more  detailed 
analysis  of  the  continuum  for   Li  in  terms  of  a  -  deuteron 
disintegration . 

The  primary  purpose  of  this  work  was  to  extend  the  de- 
tailed Saskatchewan  work  of  Ref.  5  to  the  continuum  of   Li. 
In  addition,  during  the  course  of  this  experiment  a  check 
was  made  for  a  large  inelastic  peak  observed  but  not  reported 
by  Monson  [Ref.  6] .   As  corroborating  evidence  for  this  peak 
was  not  found,  it  is  assumed  that  the  scattering  of  electrons 
off  the  target  frames  which  have  since  been  modified  was  the 
cause  of  the  apparent  peak. 


II.   THEORY 

A.   ELECTRON  SCATTERING 

When  a  high  energy  electron  is  incident  upon  a  target 
material,  it  is  scattered  by  the  Coulomb  force  developed 
between  its  negative  charge  and  the  positively  charged  nu- 
cleus.  Depending  upon  whether  or  not  energy  is  converted 
from  kinetic  to  internal  nuclear  motion,  the  scattering 
event  is  termed  elastic  or  inelastic. 

In  elastic  scattering  no  excitation  energy  is  imparted 
to  the  nucleus.   The  scattered  electron  energy  is  essentially 
the  same  as  the  incident  energy,  being  reduced  only  slightly 
by  the  recoil  of  the  relatively  massive  nucleus.   If  the  in- 
cident electron  has  a  wavelength  comparable  to  the  size  of 
the  nucleus  (i.e.  an  electron  with  energy  on  the  order  of 
100  MeV) ,  then  electron  scattering  experiments  can  be  used 
to  study  the  charge  distribution  of  the  target  nucleus.  The 
quantitative  relationship  is  given  in  terms  of  charge  form 
factors  discussed  below. 

In  inelastic  scattering  some  incident  particle  energy 
excites  the  nucleus,  as  well  as  imparting  recoil.   The  ex- 
cited nucleus  can  assume  discrete  energy  levels  and  then 
decay  to  lower  energy  levels  by  emitting  gamma  rays  or  by 
particle  emission.   Continuous  excitation  can  occur  when 
the  nucleus  absorbs  sufficient  energy  to  break  up  into 
smaller  parts  which  can  have  a  continuum  of  final  state 
energies . 
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B.   ELASTIC  SCATTERING 

The  differential  cross  section  per  unit  solid  angle  for 
a  relativistic  electron  scattering  elastically  from  a  spin- 
less  point  nucleus  of  charge  Ze  was  shown  by  Mott  [Ref.  7] 

to  be 

/dffCE  ,6)\     /da(E  ,8)\      /Ze2\2  cos2  9/2    1 


V    dft    /Mott  \    dft    /point  ^2Eo/   sin1*  9/2    R  ,  (n_i) 
where  EQ  =  incident  electron  energy 

2  F        9 

and   R  =  1  +  2.  sin    9/2   is  the  recoil  factor  for  the 

Mc2 

7       1 

target  nucleus  of  mass  M.   For  E0<<  Mc  ,  =•  - 1  and  this  term 
is  often  omitted  from  the  above  equation.   It  is  assumed 
that  the  electron  is  extremely  relativistic,  that  is,  E  >>mc2 

The  Mott  cross  section  defined  above  for  a  point  nucleus 
must  be  modified  by  a  factor  which  takes  into  account  the 
finite  charge  distribution  of  the  target  nucleus.   The  square 
of  this  charge  form  factor  is  the  ratio  of  the  observed  dif- 
ferential elastic  cross  section  to  the  Mott  cross  section. 

In  elastic  scattering  the  momentum  transfer  is  q  =  pf-pQ, 

->-      -»- 
where  pQ  and  p^  are  the  momenta  of  the  incident  and  scattered 

electrons  respectively.   It  follows  that 

2p 

■+> 


sin  9/2        2E 


1  +  ^  sin2  9/2 


x.        *c 


-in  9/2  ,  (for  E  <<Mc2 ) 


s: 


'2 


(II-2) 


and  thus  q  depends  only  upon  the  incident  electron  energy 
E   and  the  scattering  angle  9  [Ref.  8]. 

For  low  Z  elements ,  the  plane  wave  Born  approximations 
are  used  in  determining  the  cross  sections  for  the 
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excitation  of  nuclear  states.   The  plane  wave  Born  approxima- 
tion  (PWBA)  consists  of  the  following  three  approximations: 

1.  First  order  perturbation  theory  is  adequate  to  de- 
scribe the  interaction  between  the  electron  and  the  nucleus. 

2.  The  incident  and  scattered  electrons  are  treated  as 
plane  waves. 

3.  Relativistic  effects  of  the  recoiling  nucleus  are 
ignored  [Ref .  9] . 

Considering  the  incident  and  scattered  electrons  as  plane 


waves ,   ii   =  e 
o 


-ik  *  r 


■ik. 


and   ipe  =  e   f    respectively.   The 
differential  elastic  cross  section  is 


'daCE  ,6) 
o 


o  | 


da(E  ,9> 


dft 


'da  ( E  ,  9  > 
o ' 


dft 


Matt 


Mott 


Mott 


nuclear 
volume 


nuclear 
volume 


ip  -  p  Cr )  ty      d 3  r 


-ik,..r  ,~*~>.      -ik  .  r  ,  3 
e    f    p(r)e    o   dr 


_1 
Ze 


-ik -r   ,^nj3 
J        e       p (r )d  r 

nuclear 

volume 


(II-3) 


Since  k  =  k^  -  k   =  f-, 
f    o    h 


faa  CE  ,  9  )s 

o  ' 

dft 


'daCE  ,9)\ 


FCq) 


(II-4) 


/Mott 

where  F (q)  =  charge  form  factor.   When  using  the  PWBA  for 
a  spin-zero  nucleus,  the  charge  form  factor  is  thus  the 
Fourier  transform  of  the  nuclear  charge  distribution. 
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C.   INELASTIC  SCATTERING 

The  momentum  transferee!  to  the  nucleus  in  an  inelastic 
scattering  event  is  a  function  of  the  incident  and  scattered 
electron  energies  and  the  scattering  angle  (provided  that 
EQ<<  Mc2) ,  i.e., 

q2=  2  kQ2   (1  -  cos  9)  (1  -  kAQ)  +  k2 


=  2(^°- 


^P_\        (i  -  cos  6)  (l   -  f!x\   +  f^x\       ,      (II-5) 


where  kQ  =  Eo/Tic, 
k   =  E/nc, 


and   E„  =  excitation  energy. 


It  is  possible  to  derive  analytical  expressions  for  the 

differential  inelastic  electron  scattering  cross  sections 

using  the  PWBA.   With  these  approximations  the  differential 

inelastic  cross  section  (-377       can  be  written  as  a  sum 

\auj  PWBA 

over  the  cross  sections  for  electric  (e)  and  magnetic  (m) 
transitions  of  various  multipole  orders: 

/da  \  _     r*  /da\  r  (d£\ 

WpWBA    ~     W%X  X    Ka     , 

(II-6) 


where 


rdd\  2         2  A,     -2 

vdJT  j  ,    =   a   aX^     ko 


'eX 


X(X  +  l)"1B(CA,q,I    +I    )V-  C9) 


+    B(EX,q,I    ■>!    )VTC9)1  R 
^      o       x      T        J 


o      x      L 

-1 


(II-7) 


and 

/S£\  =    „a_    nzXV    -2RfMKn.T    +  T     )V_C9)R"X 


Ml  -    ^      ko      3(1!^'Io+Ix)VT(-9]R         '  CI1-81 

/r"1 


=    a    a. 
'mX 
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In  these  equations , 

(2x+d:  : 


ax  =  4tpx-1(x  +  1) 


R  =  1  +  he  (-2J      (l-cos0), 
X      =  transition  multipolarity , 
a   =  fine  structure  constant, 
EQ  =  primary  electron  energy, 
9   =  scattering  angle  and 

M  =  nuclear  mass. 

2 
If  E  > > Mc   for  the  electrons,  then 

VL(9)  =  x^(l+cos9  )/(y-cos8  )2   and  (II-9) 

VT(9)  =  ^(2y+l-cos9 )/(y-cos9 ) (l-cos9)  ,  (11-10) 

_  2 
with  y  =  1  +  ^x  /2EQ  (EQ  -  Ex) . 

For  E  >>  E„  we  have  y  -*-l  resulting  in 

VL(9)  =  cos-2  9/2  /  usin*  9/2  (11-11) 

and 

VT(9)    =    1+sin2  9/2    /    8  sin  *  9/2 

=    cos2    I   (^+tan2    |o    /    Usin"    |  (11-12) 

The  R  factor  accounts  for  the  effect  of  the  recoiling 
target  nucleus  on  the  final  density  of  states.   The  func- 
tions VL (  9  )  and  VT (  9 )  refer  to  the  longitudinal  and  trans- 
verse components  of  the  four-vector  nuclear  current  density 
resolved  into  components  parallel  and  perpendicular  to  the 
three-vector  momentum  transfer  q.   The  longitudinal  (paral- 
lei  to  q)  component  represents  an  interaction  of  the  elec- 
tron with  the  nuclear  charge.   The  transverse  (perpendicular 
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to  q)  component  represents  the  contributions  of  the  nuclear 
current  and  magnetization  densities.  B(  X  ,  q,I  -»-I  )  repre- 
sent the  reduced  nuclear  transition  probabilities  [Ref .  10] . 

The  different  angular  functions  VL (  © )  and  VT (  0 )  allow 
one  to  emphasize  one  type  of  transition  over  another.   At 
backward  (9>  90°)  scattering  angles  V_ (  8  ) >  VT  (  9  ) ,  while 
at  forward  scattering  angles  V  (  9  )  predominates.   Thus,  mag- 
netic transitions,  which  only  involve  V  (  9  )  are  readily 
identified  by  scattering  at  180°,  while  electric  multipole 
transitions  dominate  at  forward  scattering  angles. 

For  low  Z  elements  the  PWBA  given  sufficiently  accurate 
estimates  of  the  inelastic  differential  cross  sections.   Dis- 
torted wave  Born  approximations  (DWBA)  therefore  will  not  be 
discussed  in  this  paper. 

In  First  Born  approximation,  equations  (11-11  and  12)  re- 
duce equations  (II-7  and  8)  to 


'd8\  /d  9\ 


1  I 


/Mott 
and      {d£\        =    fd£]  r_  Q 


M        VdVMott   L(^  +  tan  2)TT    ^] 


(11-13 


(11-14) 


respectively,  where  F   (q)  and  F  ^(q)  are  the  inelastic  form 
factors  for  the  longitudinal  and  transverse  photon  absorption 
[Ref.  11].   By  maintaining  a  constant  q  for  different  scatter- 
ing angles  by  varying  the  incident  energy,  it  is  possible  to 

2  2 

determine  F   (q)  and  F   (q)  separately. 

D.   INELASTIC  FORM  FACTORS 

In  order  to  determine  the  inelastic  form  factors  one 
must  start  by  calculating  the  area  under  the  elastic  peak. 
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Two  types  of  background  must  be  removed  from  the  inelastic 
spectrum.   These  include  electrons  which  are  originally- 
scattered  by  the  target  and  which  subsequently  scatter  off 
the  spectrometer  walls  and  enter  the  detector  system  and  the 
elastic  peak  radiation  tail  caused  by  photon  emission  in 
conjunction  with  nuclear  scattering,  plus  the  effects  of 
ionization  and  straggling. 

If  one  assumes  that  only  one  photon  is  emitted  in  conjunc- 
tion with  the  nuclear  scattering  event  and  that  the  PWBA  is 
applicable,  the  contribution  to  the  radiation  tail  caused  by 
photon  emission  can  be  reliably  evaluated.   The  determination 
of  the  effects  of  stragling  are  less  precise  and  one  must  re- 
normalize  the  entire  radiation  tail  to  fit  the  experimental 
spectrum  in  the  region  of  interest  [Ref .  5] . 

After  extending  the  shape  of  the  elastic  radiation  tail 
into  the  inelastic  region  and  subtracting  it  from  the  ex- 
perimental spectrum,  the  inelastic  spectrum  remains.   The 
differential  inelastic  cross  section  of  a  narrow  inelastic 
peak  can  be  determined  by 

(do  (E  ,  9  )\         _      in  ( o  '   \ 


\       /in     el  \        / 

=  Ain  /da(Eo,9)\ 

"  *e7  \        ^    Aott  '  F2(q°}  '      (H-15) 

/dcrCE  ,9)\ 
where  (  2 )        =  differential  elastic  cross  section, 

\      dQ,  J 

A.  =  area  of  inelastic  peak, 

in  c 

A  ,  =  area  of  elastic  peak, 


lel 

'da  (.: 


=  differential  elastic  Mott  cross 


d^   /       section  and 
/Mott 
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2 

F  Cq  )        3  elastic  form  factor 
■o 

2 
The  inelastic  form  factor  F  Cq)  can  then  be  determined  by 


da  CE  ,  9  )) 
FZ(q)  =  V  dfl  J  (11-16) 

/da(E  ,8X 


\       Mott 

Equations  (11-15)  and  (11-16)  assume  a  common  shape  be- 
tween the  elastic  and  inelastic  peaks  and  that  the  radiation 
corrections  are  the  same  for  both  elastic  and  inelastic 
scattering.   The  inelastic  cross  section  is  expressed  as  an 
area  per  unit  solid  angle  since  it  is  determined  in  terms  of 
an  integral  over  the  inelastic  peak. 

In  determining  the  inelastic  cross  section  for  the  exci- 
tation of  a  continuum  of  final  states  the  cross  section  is 
expressed  in  terms  of  an  area  per  unit  final  energy  per  unit 
solid  angle.   After  determining  the  area  (A.  )  by  integrating 
over  the  desired  energy  range,  the  continuum  cross  section  can 
be  determined  by  equation  (11-15).   The  continuum  form  factor 


can  then  be  determined  using  equation  (11-16) 


Ref .  3 
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III.   EXPERIMENTAL  EQUIPMENT  AND  PROCEDURE 

A.   NAVAL  POSTGRADUATE  SCHOOL  LINEAR  ACCELERATOR 

The  NPS  120  MeV  LINAC  is  an  electron  accelerator  similar 
to  the  Stanford  University  Mark  III  linear  accelerator  de- 
scribed by  Chodorow  et  al  [Ref .  13] .   The  NPS  LINAC  consists 
of  three  10  foot  accelerator  sections  and  is  powered  by 
three  klystrons,  each  capable  of  22  megawatts  peak  power. 
The  machine  is  pulsed  at  a  rate  of  60  pulses  per  second  with 
each  beam  pulse  having  approximately  one  microsecond  dura- 
tion.  The  maximum  average  beam  current  is  approximately  one 
microampere. 

Beam  deflection  magnets  and  energy  defining  slits  are 
used  to  select  the  beam  energy,  which  can  range  from  3  0  to 
105  MeV  with  energy  resolution  of  0.3%.   The  energy  resolu- 
tion attained  in  the  experiments  reported  herein  was  0.6%. 
Target  thickness  is  the  main  reason  for  the  lower  resolution 
For  this  reason,  the  478  KeV  inelastic  peak  could  not  be 
resolved;  therefore,  its  area  is  added  to  that  of  the  elas- 
tic peak. 

The  field  of  the  beam  deflection  magnet  could  be  meas- 
ured by  a  nuclear  magnetic  resonance  (NMR)  probe,  but  this 
probe  did  not  function  during  these  runs.   A  total  of  four 
quadrupole  magnets  are  available  to  focus  the  beam  to  a 
diameter  of  approximately  two  millimeters  at  the  target. 
In  recent  experiments  only  the  last  two  have  been  used, 
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because  the  first  pair  cause  the  beam  to  expand,  hit  the  beam 
pipes  and  cause  excess  background. 

B.   SPECTROMETER  AND  COUNTING  SYSTEM 

The  NPS  LINAC  has  a  movable,  16-inch  radius,  double  focus- 
ing magnetic  spectrometer  capable  of  measuring  electrons 
scattered  at  angles  from  45  to  150  degrees  in  15  degree  steps. 

The  energy  distribution  of  scattered  electrons  is  measured 
by  varying  the  field  of  the  spectrometer  magnet.   Ten  NE102 
plastic  scintillation  counters  located  in  the  focal  plane  of 
the  spectrometer  serve  as  detectors  for  the  ten  counting 
channels.   The  spread  of  these  counters  span  a  range  ±   1.5 
percent  about  the  spectrometer  energy.   Each  counter  thus 
detects  an  energy  spread  of  0.3  percent  of  the  central  spec- 
trometer energy.   Two  large  detectors,  each  spanning  the 
same  field  width  as  the  ten  smaller  detectors,  act  as  dual 
coincidence  backing  counters.  An  electron  passing  through  a 
front  detector  and  both  backing  detectors  will  record  one 
count  on  the  Jorway  Model  183  6  Dual  Channel  Scalar.   The 
accelerator  gun  pulse  and  the  counter  gating  circuit  are 
synchronized  so  that  counts  are  recorded  only  during  pulses. 
This  prevents  background  counts  from  being  counted  between 
pulses . 

The  scalars  are  capable  of  recording  periodic  pulses  up 
to  150  MHZ.   However,  since  the  beam  is  pulsed,  the  count 
rate  was  limited  to  50  counts  per  second  in  the  counter 
registering  the  elastic  peak  and  150  counts  per  second  in 
the  back  counters  while  scanning  the  inelastic  area  of  the 
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spectrum.  These  count  rates  avoided  significant  count  rate 
losses.   A  Faraday  cup  is  capable  of  calibrating  the  charge 
collected  on  the  Secondary  Emission  Monitor  (SEM) .   Past 
calibrations  showed  less  than  0.1%  change  in  the  SEM  cali- 
brations . 

C .   TARGETS 

7 
Li  targets  with  thicknesses  of  106.0,  147.5,  and  150.0 

mils  were  used.   Since  lithium  reacts  quickly  with  air  to 
form  lithium  nitride,  target  preparation  had  to  be  conducted 
with  the  lithium  immersed  in  an  oil  bath.   The  targets  were 
hydraulically  pressed  between  polished  stainless  steel 
plates  separated  by  metal  spacers  to  approximate  the  desired 
thicknesses.   This  method  precludes  preparing  targets  ac- 
curately to  a  predetermined  thickness.   Several  micrometer 
readings  of  the  target  thickness  were  taken  over  the  surface 
of  the  target.   Thickness  readings  varied  by  less  than  0.3 
percent.   The  targets  were  mounted  in  an  aluminum  ladder 
and  the  target  chamber  was  quickly  evacuated  to  approximate- 
ly 1.5  x  10    torr  to  limit  target  contamination.   Chamber 
vacuum  had  to  be  released  when  changing  the  scattering  angle 
and  the  targets  thus  became  exposed  to  air  for  short  periods 
of  time.   The  effect  of  this  is  seen  in  Figures  IV-4,  6,  7 
and  8  which  show  peaks  on  the  high  energy  side  slope  of  the 

elastic  scattering  peak.   This  was  caused  by  the  difference 

7 
in  recoil  of  the   Li  nuclei  and  the  N  nuclei.   Elastic  peak 

areas  were  adjusted  to  account  for  this  effect. 
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D.   DATA  ACCUMULATION 

Beam  energy  and  scattering  angle  were  kept  constant  dur- 
ing each  data  accumulation  run.   The  spectrometer  energy  was 
decreased  in  50  KeV  steps  over  the  elastic  peak  and  in 
0.1  MeV  steps  over  the  inelastic  region  of  interest.   The 
integration,  capacitance  and  beam  intensity  were  chosen  to 
prevent  exceeding  maximum  count  rates  discussed  in  Section 
III-B.   For  the  90  MeV  runs  the  spectrometer  energy  was 
initially  set  at  100  MeV  and  then  reduced  to  a  starting  point 
approximately  two  MeV  above  the  elastic  peak.   Elastic  data 
was  recorded  in  50  KeV  steps  to  approximately  six  MeV  below 
the  elastic  peak.   During  the  inelastic  runs  the  beam  was 
retuned  to  give  maximum  intensity  consistent  with  the  above 
mentioned  count  rate  restrictions  and  the  spectrometer  energy 
was  returned  to  a  starting  position  which  showed  the  elastic 
peak  maximum  registering  in  counting  channel  number  ten.   In 
order  to  follow  the  same  hysteresis  loop,  the  spectrometer 
magnet  energy  was  returned  to  100  MeV  and  then  reduced  to 
the  starting  point  energy.   Inelastic  data  was  taken  in  0.1 
MeV  steps  down  to  a  point  approximately  nine  MeV  below  the 
elastic  peak.   A  similar  procedure  was  used  for  the  105  MeV 
runs.   In  this  case  the  spectrometer  energy  was  initially 
set  at  110  MeV  and  then  reduced  to  a  starting  point  approxi- 
mately 2.5  MeV  above  the  elastic  peak.   Elastic  data  was  re- 
corded in  50  KeV  steps  to  approximately  seven  MeV  below  the 
elastic  peak.   Prior  to  collecting  the  inelastic  data  the 
spectrometer  energy  was  returned  to  110  MeV  and  then  to  the 
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inelastic  starting  point  as  per  the  arguments  presented 
above.   Inelastic  data  was  recorded  in  0 . 1  MeV  steps  down 
to  a  point  approximately  ten  MeV  below  the  elastic  peak. 

Incident  energy,  integration,  machine  energy,  capaci- 
tance, elapsed  time  and  the  counts  for  the  ten  counting  chan- 
nels and  two  back  counters  were  automatically  recorded  on  a 
digital  readout.   This  data  was  electrically  transmitted  to 
an  IMSAI  8080  micro-computer  and  displayed  on  an  ancillary 
cathode  ray  tube  (CRT)  where  it  could  be  verified  before 
being  recorded  on  both  cassette  tape  and  magnetic  seven- 
track  tape.   A  teletype  also  provided  a  hard  copy  of  re- 
corded data.   The  magnetic  tape  was  used  to  provide  input 
to  Code  ASC  described  in  Appendix  A.   The  hard  copy  was  used 
in  editing  the  data  after  it  had  been  transferred  from  the 
magnetic  tape  to  disk  storage  in  the  NPS  IBM  36  0  computer. 
The  cassette  tape  provided  a  backup  and  could  be  used  to 
re-record  a  seven- track  tape  if  necessary. 

Code  ASC,  which  uses  the  input  from  the  spectrometer 
energy,  counting  channels,  integration,  capacitance  and 
elapsed  time  of  each  run  to  compute  a  single  spectrum, 
automatically  corrects  the  energy  of  each  counter  to  cor- 
respond to  the  true  energy  being  detected.   Relative  count- 
er efficiencies  are  therefore  taken  into  account  because 
each  point  in  the  final  spectrum  is  actually  the  average 
of  noints  from  each  counter. 
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IV.   DATA  ANALYSIS 

A.   ELASTIC  SPECTRUM 

The  spectrum  generated  by  Code  ASC  was  plotted  and  the 
area  under  the  elastic  peak  was>  determined  as  described  be- 
low.  Integration  of  the  area  under  the  elastic  peak  was  car- 
ried out  from  the  upper  limit  of  the  peak  to  an  energy  two 
full  widths  at  half  maximum  (FWHM)  down  from  the  peak  energy 

The  area  of  the  elastic  peak  thus  obtained  must  be  cor- 
rected for  virtual  and  real  emission  and  absorption  of  radi- 
ation during  scattering.   Schwinger  [Ref.  14]  originally 
developed  a  multiplicative  correction  for  virtual  emission 
and  absorption  assuming  a  recoiless  target  nucleus.   As 
modified  by  Tsai  [Ref.  15]  to  include  recoil,  thus  multi- 
plicative constant  is 


K   =  e 
s 


(IV-1) 


where  5   -  2a 

s 

7T 


In   1 


13 

72    12 


2E.Mc 
1 


—  (i4i) 


+  17 

35"  ,  (IV-2) 


AER 

V 
where   a  =  the  fine  structure  constant, 

E.   =  the  incident  electron  energy, 

E  =  elastic  peak  energy  (E  )  -  lower  cutoff  energy 
used  in  determining  elastic  peak  area, 

R   =  recoil  factor  previously  defined, 

M   =  rest  mass  of  target  nucleus  and 

m   =  electron  rest  mass. 


23 


Bethe  and  Ashkin  [Ref.  16]  developed  the  correction  for 
the  emission  of  real  radiation 

Kb  =  e  D  ,  (IV-3) 

where  *b  =  !o   ln  fi.      ,  (IV-4) 

ln2    AER3/2 
t   =  the  effective  target  thickness  in  radiation 
lengths , 

and  the  other  parameters  are  as  defined  for  equation  IV-2 

above.   Equation  IV- 4  is  valid  only  if 

The  elastic  area  then  becomes 

A(V  ■  KsKb  A(V  observed.  (IV"5> 

Although  background  radiation  was  subtracted,  the  area 

due  to  the  unresolved  478  KeV  excited  state  was  not;  thus 
the  elastic  area  are  slightly  larger.   The  differential 
elastic  scattering  cross  section  was  determined  using  equa- 
tion II-4.   In  this  equation  the  F(q)  used  was  as  reported 
by  Bumiller  et  al  [Ref.  17]  as 

F(q)  =  d(l-cq  )e"aq  ,  (IV-2) 

where  d   =1.0, 

c   =0.16  9  and 

a   =  0.713. 

B.   CALCULATION  OF  THE  ELASTIC  PEAK  RADIATION  TAIL 

The  radiation  tail  is  primarily  the  result  of  radiation 
during  scattering  and  radiation  before  or  after  scattering 
in  the  electromagnetic  field  of  a  different  nucleus.   The 
effect  of  radiation  during  scattering  can  be  evaluated,  as 
shown  by  Barber  et  al  [Ref.  12] ,  using  the  Schiff  [Ref.  18] 
calculation  of  the  differential  cross  section, 
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d*.aa(E)  fd0(Eo)^         fx  +  e£\      f"    ln  /2E    sin9/2\    _    l" 


dfi  \     dft       /  \  o7 

1    +    ^O2         F2(q) 


(]VE>        V  T2  F2(qo)y 


(IV-7) 


where     /da(EQ)\     =  /da(EQ)\  .      F2    (qQ)     , 


\   dfl      y      \   da    /  Mott 

a  =  fine  structure  constant  and 

m  =  electron  mass. 
The  area  of  the  elastic  peak  A(E  )  was  determined  in 
counts  -  MeV  units.   The  spectrum  of  the  Schiff  radiation 
tail  N  (E)  in  counts  was  determined  by 

(  da_    (E0)\ 
A(EQ)     V  dft      J 


NS(E)     /  d^  aa(E) 


2  „  /^N    or 


dfl 


_   A(E0) 


d^  aa(E)N 


NS(E) 

'dg     (E0T 

dfi 


(IV-8) 


E 
where  the  factor  ■=—  corrects  for  counter  resolution. 

Eo 

The  effect  of  radiation  before  or  after  scattering  due 
to  brehmsstrahlung  is 

d2  crB(E)   =  bl      1     fdg  (EQ)     da  (E)_ 
dH  dE        2     (Eo~E)    [   dQ     +  ~7^   J  ,  (IV-9) 

where  b     (target  thickness  in  radiation  lengths) 
1  ln2 

• 

The  spectrum  of  the  bremsstrahlung  radiation  tail  N,  (E) 

was  then  computed  in  a  manner  similar  to  equation  IV-8  above 

and  the  sum  of  N_  (E)  and  N'D(E)  yield  the  spectrum  of  the 

total  radiation  tail. 
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The  radiation  tail  then  had  to  be  normalized  as  discussed 
earlier  to  fit  the  area  of  interest  on  the  elastic  peak's 
lower  energy  slope.   The  tail  was  normalized  from  1  MeV  ex- 
citation energy  to  the  a  -triton  threshold  at  ^  2.4  5  MeV 
excitation  energy.   The  normalized  radiation  tails  are  de- 
picted in  Figures  IV- 1  through  8.   The  noramlized  radiation 
tail  was  then  subtracted  from  the  raw  data  spectrum  to  ob- 
tain the  inelastic  continuum  as  shown  in  Figures  IV- 9  through 
16. 

C.   INELASTIC  SPECTRUM 

In  order  to  separate  the  inelastic  peak  at  4.55  MeV  from 
the  continuous  inelastic  spectrum,  a  fitting  routine  was 
used  to  generate  a  line  fit  to  the  data  on  either  side  of 
the  inelastic  peak  from  threshold  to  6  MeV  excitation  energy. 
This  line  fit  is  represented  by  the  dashed  line  in  Figures 
IV- 9  through  16.   The  line  fit  was  subtracted  from  the  spec- 
trum in  order  to  yield  the  contribution  of  the  inelastic 
peak.   The  area  of  the  inelastic  peak  was  then  determined 
and  compared  to  the  corrected  elastic  peak  area  as  shown  in 
Table  IV- 2.   Per  arguments  presented  in  Section  II-D,  the 
elastic  and  inelastic  peaks  were  assumed  to  have  the  same 
shape  and  thus  the  inelastic  cross  section  can  be  determined 
by  equation  (11-15).   The  results  are  discussed  in  Section 
V-B  and  shown  in  Table  IV- 2. 

The  inelastic  continuum  cross  section  per  unit  final 
energy  was  determined  by  integrating  the  fitted  line  shape 
described  above  from  the  threshold  excitation  energy  as  a 
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lower  limit  to  the  desired  upper  limit  <_  6  MeV  excitation 
energy.   The  corresponding  inelastic  form  factors  for  both 
the  continuum  and  the  inelastic  peak  were  determined  by 
dividing  the  computed  integrated  cross  sectons  by  the  Mott 
cross  sectons  evaluated  at  the  appropriate  incident  ener- 
gies.  Results  are  presented  in  Tables  IV- 3  and  4. 

D.   ERROR  ANALYSIS 

Since  the  experimental  cross  sections  and  form  factors 
were  determined  by  the  ratio  of  elastic  and  inelastic  peak/ 
continuum  areas  resulting  from  data  accumulated  during  a 
single  run,  many  of  the  systematic  errors  cancelled;  e.g. 
target  thickness,  solid  angle,  SEM  efficiency  and  spectro- 
meter efficiency. 

As  mentioned  earlier,  the  spectra  of  certain  runs  showed 
evidence  of  nitrogen  contamination  which  added  to  the  elas- 
tic peak  area.   The  uncertainty  in  adjusting  the  elastic 
peak  areas  due  to  this  effect  was  estimated  to  be  less  than 
0.2%. 

The  inability  to  resolve  the  47  8  KeV  level  in  this  ex- 
periment means  that  the  elastic  peak  areas  are  too  large. 
Using  the  arguments  presented  by  Suelze  et  al  [Ref.  19]  it 
is  estimated  that  the  error  introduced  will  be  approximately 
two  percent  at  most. 

The  statistical  error  computed  by  Code  ASC  and  depicted 
in  Figures  IV- 1  through  8  show  that  these  errors  were  in- 
significant in  determining  the  elastic  to  inelastic  peak 
area  ratios.   Relatively  small  errors  would  also  be  introduced 


27 


by  the  step  sizes  used  in  the  experiment  and  the  limitations 
which  they  impose  on  locating  the  elastic  and  inelastic  peak 
energies.   This  will  have  a  correspondingly  small  effect  on 
determining  the  spectrum  of  the  radiation  tail  of  the  elastic 
peak. 

The  largest  error  is  introduced  by  the  combined  effect 
of  normalizing  the  radiation  tail  to  the  areas  of  the  elas- 
tic peak  as  already  described  in  Section  IV-B,  and  in  gen- 
erating a  least  squares  fit  to  the  inelastic  continuum  from 
threshold  to  6  MeV  excitation  energy.   The  worst  case  would 
occur  at  the  highest  q  value  since  the  integrated  cross  sec- 
tion of  the  continuum  decreases  as  q  increases.   Since  the 
line  shape  spectrum  of  the  fitted  curve  is  subtracted  from 
the  inelastic  peak,  errors  which  tend  to  increase  the  area 
of  the  continuum  have  the  reverse  effect  on  the  area  of  the 
inelastic  peak  and  vice  versa.   It  is  estimated  that  these 
effects  could  introduce  errors  into  the  determined  inelas- 
tic cross  sections  and  form  factors  of  as  much  as  -5%. 
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TABLE  IV- 1 
ELASTIC  SCATTERING  PARAMETERS  FOR  7Li 


q\ 

(d^Mott 

\anj 

(MeV) 

e 

(F-2) 

F(q) 

(mb/sr) 

(mb/sr) 

89.47 

60 

.2000 

.8378 

3.871X10"1 

2.717xl0-1 

105 

60 

.2682 

.7885 

2.881X10"1 

1.791xl0_1 

89.22 

75 

.2935 

.7709 

1.484xl0_1 

8.819xl0"2 

89.22 

90 

.3902 

.7072 

6.527xl0"2 

3.264xl0"2 

105 

75 

.3969 

.7030 

1.094xl0_1 

5.406xl0~2 

89.22 

105 

.4868 

.6486 

3.060xl0"2 

1.287xl0"2 

105 

90 

.5251 

.6267 

4.830xl0"2 

1.897xl0"2 

105 

105 

.6501 

.5600 

2.280xl0"2 

7.147xl0"3 
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TABLE  IV-2 
INELASTIC  PEAK  SCATTERING  PARAMETERS  FOR  ?Li 


fda\  fda\ 

2  Uo/M  **  vein ) 

So  a   /A  V  /Mott  V  / 

(F  l)                in   el  (mb/sr)  F(q)          (mb/sr) 

.1898  2.374xl0~3  3.371X10"1  4.032xl0"2  6.451xl0~4 

.2565  1.065xl0"2  2.881xl0_1  3.138xl0"2  1.903xl0~3 

.2735  5.986xl0"3  1.434x1a"1  5.965xl0"2  5.279x1c"11 

.3703  1.104xl0"2  5.527xl0"2  7.224xl0"2  3.406xl0"4 

.3793  3.377xl0"3  1.094X10"1  6.434xlQ"2  4.523xl0"4 

.4522  1.757xl0"2  3.060xl0"2  8.593xl0~2  2.262xl0~4 

.5018  1.918xl0"2  4.330xl0~2  3.679xl0'2  3,638xl0"4 

.6227  2.871xl0"2  2.280xl0"2  9.487xl0"2  2.052xl0~4 
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TABLE  17-3 

CONTINUUM  PARAMETERS 
(FROM  THRESHOLD  TO  INELASTIC  PEAK  ENERGY) 


2 

^2 
CF  Z) 

Integrated 

A    .  /A  , 
cont   el 

4.  308xl0"3 

Cross  Section 

Cmb/sr ) 

F(q) 

.  1898 

1.171xl0"3 

5.499xl0"2 

.2565 

4.239xl0"3 

7.592xl0"4 

5 .134xl0"2 

.2735 

3.565xl0"3 

7.553xl0"4 

7.135xl0"2 

.  3703 

1.091xl0"2 

3. 562xl0"4 

7. 338xl0"2 

.  3793 

1.069xl0~2 

5. 778X10"4 

7.268xl0"2 

.4622 

1. 171xl0"2 

1.  507xl0-14 

7.013xl0"2 

.5018 

1.548xl0"2 

2 .936xl0"4 

7. 797xl0"2 

.6227 

1. 571xl0"2 

1.123xl0"4 

7.017xl0~2 
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TABLE  IV  -4 

CONTINUUM  PARAMETERS 
(.FROM  THRESHOLD  TO  6  MeV  EXCITATION  ENERGY 


2 

(p32> 

Integrated 

A    . /A  . 

cont   el 

1.075xlQ~2 

Cross  Section 

Cmb/sr ) 

F(q) 

.  1865 

2.922xl0"3 

8.633xl0"2 

.2526 

1. 175xl0"2 

2 .104xl0~3 

8.545xl0"2 

.2738 

2.013xl0"2 

1. 775xl0"3 

1.094X10"1 

.  3638 

2.644xl0"2 

8.630xl0"4 

1.150xl0-1 

.  3739 

2. 380xl0"2 

l,237xl0"3 

1.085X10"1 

.4539 

3.175xl0~2 

4.087X10"14 

1.156X10'1 

.4946 

3.744xl0~2 

7,102xl0"4 

1. 213x1a"1 

.6122 

4.619xl0"2 

3,  301xl0'4 

1.203xl0_1 
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V.   RESULTS  AND  CONCLUSIONS 

A.  COLLECTIVE  RESULTS 

Results  of  eight  separate  data  runs  (four  at  a  primary 
electron  energy  of  ^  90  MeV  and  four  at  a  primary  electron 
energy  of  105  MeV  and  at  scattering  angles  of  60,  75,  90  and 
105  degrees)  are  depicted  in  Figures  IV- 1  through  19  and 
Tables  IV- 1  through  4.   Due  to  an  experimental  energy  reso- 
lution of  0.6%,  the  inelastic  peak  at  478  KeV  could  not  be 
resolved.   The  figures  show  an  inelastic  continuum  starting 
in  the  vicinity  of  the  a  -triton  threshold  (2.45  MeV  excita- 
tion energy)  and  extending  to  higher  excitation  energies.  An 
elastic  peak  which  has  previously  been  reported  at  4.6  3  MeV 
[Ref.  2,    3,  4,,  19]  was  found  at  4.55  MeV  excitation  energy 
and  its  spectrum  is  superimposed  on  the  inelastic  continuum. 

The  curves  for  inelastic  peaks  and  continuum  cross  sec- 
tions are  plotted  separately  for  the  different  energies  be- 
cause form  factors  are  not  functions  of  q  only  if  transverse 
multipoles  are  involved. 

B.  4.55  MeV  INELASTIC  PEAK 

As  mentioned  above,  the  4.55  MeV  inelastic  peak  (4.63 
MeV)  has  been  previously  reported  and  analyzed  on  several 
occasions.   Within  the  primary  electron  energy  groups  dis- 
cussed above  the  ratio  of  the  inelastic  peak  area  to  the 
elastic  peak  area  and  the  inelastic  form  factors  increase 
as  q  increases.   The  differential  Mott  and  inelastic  cross 
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section  decrease  as  q  increases.   Experimental  values  are 
presented  in  Table  IV- 2  and  plotted  in  Figure  IV-18.   The 
data  run  at  105  MeV  with  a  scattering  angle  of  60°  is  an 
exception  to  these  arguments  and  it  will  be  discussed  sep- 
arately below. 

C.  INELASTIC  CONTINUUM 

As  shown  in  Tables  IV- 3  and  4,  the  ratio  of  the  continu- 
um areas,  integrated  from  threshold  to  4.55  MeV  excitation 
in  the  case  of  Table  IV- 3  and  to  6  MeV  excitation  in  the 
case  of  Table  IV- 4,  to  the  elastic  peak  area  increases  as  q 
increases  in  each  primary  electron  energy  group.   The  differ- 
ential cross  sections  integrated  over  these  limits  decrease 
as  q  increases.   As  shov/n  in  Figure  IV-19,  the  continuum  form 
factors  increase  with  q  to  a  maximum  and  then  decrease  with 
increased  q. 

D.  DISCREPANCY 

As  depicted  in  Table  IV-2  and  Figure  IV-18,  the  ratio  of 
the  inelastic  peak  area  to  the  elastic  peak  area  for  the 
105  MeV,  60°  run  does  not  agree  with  values  expected  from  an 
analysis  of  the  other  105  MeV  runs.   In  spite  of  this,  the 
continuum  cross  section  and  form  factors  agree  closely  to 
those  values  predicted  by  extending  the  least  squares  line 
fit  of  the  other  105 MeV  runs  down  to  the  lower  momentum 
transfer  of  this  data  run.   The  cause  of  this  discrepancy 
has  not  been  determined  and  should  be  investigated  in  the 
future.   The  parameters  from  this  data  point  were  not  used 
in  generating  the  line  fits  of  Figures  IV-18  and  19. 

53 


E.   CONCLUSIONS 

Results  of  these  experiments  are  in  general  agreement 
with  results  achieved  by  Bishop  and  Bernheim  [Ref.  3]  for  a 

different  range  of  momentum  transfer.   The  continuum  observed 

.  .  7 

appears  to  be  associated  with  the  disintegration  of   Li  into 

an  alpha  particle  and  a  triton.   Data  presented  herein  should 

be  further  analyzed  in  an  effort  to  draw  conclusions  in  terms 

of  various  nuclear  models  and  multipole  transitions. 
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APPENDIX  A 
COMPUTER  PROGRAMS 

A.   ASC 

Code  ASC  is  a  revision  of  the  Code  KIK  described  in  de- 
tail by  Warshawsky  and  Webber  [Ref.  20].   This  revision  al- 
lows up  to  599  data  sets  to  be  processed.   The  input  para- 
meters of  each  data  set  are  discussed  below.   ASC  compiles 
the  separate  spectra  of  the  ten  counting  channels  into  one 
combined  spectrum  showing  energy  (MeV) ,  count  rate  (counts/ 
microcoulomb)  and  error  (counts/microcoulomb) . 

The  set  up  of  a  FORTRAN  card  deck  to  use  Code  ASC  con- 
sists of  the  appropriate  control  cards,  four  parameter  cards, 
the  data  deck  and  two  trail  cards.   When  using  the  Control 
Program/Cambridge  Monitor  System  (CP/CMS)  at  the  NPS  Computer 
Facility,  two  files  are  required.   A  FT04F001  file  serves  as 
the  data  set  file  and  a  FT05F001  file  serves  as  the  para- 
meter file.   In  the  description  which  follows  the  numbered 
headings  refer  to  separate  parameter  cards  in  the  FORTRAN 
deck  or  separate  lines  in  FT05F001  as  appropriate.   FORTRAN 
formats  appear  in  parentheses.   Parameters  marked  by  an  as- 
terisk are  optional  and  may  be  left  blank. 

1.  TEXT      Identifying  alphanumeric  remarks       (2  0A4) 

2.  ENGUP,  ENGLOW,  DENG,  TAU ,  SUMCHK  (5F10.5) 

ENGUP     Upper  spectrometer  energy  limit  for  the 
spectrum  in  MeV. 

ENGLOW    Lower  spectrometer  energy  limit  for  the 
spectrum  in  MeV. 
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DENG 

TAU* 
SUMCHK* 


Spectrometer  energy  step  width  in  MeV. 
ASC  interpolates  to  generate  missing 
data  sets. 

Counting  system  dead  time/pulse  length. 

=  0.   ASC  will  then  inform  user  when  the 
total  counter  sum  is  incorrect  in  the 
particular  data  set. 


N(l) 


,N(10)  (1014) 

Spectrometer  counting  system  normalization 
factors  relative  to  counter  number  6 .  The 
N(i) 's  are  defined  as 


N(i)  = 


Energy (6)  _  , 
Energy (i) 


x  10 


LU,  LO,  INV,  FKXAD,  FKEN 


(3I10,2F10.5 


LU*,  LO*  The  spectrum  of  counter  J,  LU   J   LO  will 
be  punched  out  onto  cards. 

INV  If  INV  0,  then  the  order  of  the  input  data 
sets  is  reversed;  i.e.  from  lower  to  higher 
energies . 

FKXAD    When   0  the  input  capacitor  voltage  (at 

full  scale)  will  be  scaled  to  units  of  10""  3 
volts  with  multiplication  by  this  factor. 
For  the  NPS  LINAC,  FKXAD  =  0  defaults  and 
the  voltage  is  recorded  in  millivolts. 

FKEN     The  input  energy  will  be  scaled  to  units 

of  KeV  when  multiplied  by  this  factor.  For 
the  NPS  LINAC,  FKEN  =  1.014. 

Data  cards  now  follow  with  two  cards  comprising  a  data  set. 

Input  parameters  are  in  the  order  recorded  on  the  magnetic 

tape  as  described  in  Section  II-D  and  are  as  follows: 

Card  1:   XHZ , MIVOLT , EO , MUEF ,XIP , XAD,DUM1, ZESS (1) . . . 
ZESS(3)  (10F7.0) 

Card  2:   ZESS ( 4) . . . ZESS ( 10) , UNTZ , DUM2 , SUM10 (10F7 . 0) 

XHZ      Spectrometer  energy  setting  (MeV) 

MIVOLT   Full  scale  voltage  setting  of 
capacitor  (mv) 
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EO       Incident  accelerator  energy  (MeV) , 

MUEF     Capacitance  of  collecting  capaci- 
tor (I-lF)  . 

XIP      Elapsed  time  (sec) , 

XAD      Final  voltage  on  the  cpacitor 
(mv)  , 

DUMl     Spectrometer  "accidentals"  counter, 

ZESS(I)   Counts  registered  in  "ith"  counter, 

UNTZ      Spectrometer  backing  counter  number 

11  reading, 

DUM2      Spectrometer  backing  counter  number 

12  reading. 

SUM10     Total  counts  registered  in  counters 
1  through  10. 

5.   Trail  cards  (10F7.0) 

(Same  format  and  parameters  as  data  cards  in  4  above) . 

When  a  card  is  encountered  with  columns  1-7  blank,  it  and  the 

next  card  are  considered  trail  cards  and  parameter  values  for 

ZESS(I)  are  given  the  following  meaning: 

ZESS(I)   If  j&   0,  then  counter  number  will  not  be  con- 
sidered when  compiling  the  total  spectrum. 
The  value  of  N(I)  on  parameter  card  3  above 
must  also  =  0. 

When  using  CP/CMS  a  FT04F001  file  is  used  as  a  data  file 
and  the  trail  cards  are  entered  as  two  additional  lines  at 
the  end  of  this  file.   All  formats  are  the  same. 

A  sample  of  the  parameter  cards  and  data  deck  for  a  hy- 
pothetical run  at  90  MeV  machine  energy  with  a  scattering 
angle  of  75  degrees  follows.   The  dashed  line  separates  the 
parameter  cards  (FT05F001  file)  from  the  data  deck  (FT04F001 
file)  for  purposes  of  clarification  only.   The  two  blank 
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trail  cards  at  the  end  of  the  data  deck  indicate  that  all 
counters  were  used. 

SAMPLE 
LI  90  EL  (90  MEV,  100  SLITS,  75  DEC,  14  4 . 0MG/CM**2)  6  JAN  79 
90.0      81.5        .05 

2050  1600  1130   730    350    0  -330  -730-1020-1390 

1.014 

090.00  001000  090.00  0054.3  030122  000000  000007  000004  000004 
000003  000004  000002  000001  000000  000000  000782  000542  000025 
. . .  (Rest  of  data  deck  follows)  . . . 
(Blank  trail  card) 
(Blank  trail  card) 

B.   HP  98  45A  PROGRAMS 

Simple  computer  programs  were  written  for  use  on  a 
Hewlett-Packard  Model  9845A  Computer.   These  programs  stored 
the  spectrum  output  from  Code  ASC  onto  cassette  tape  data 
files  with  separate  data  files  for  the  spectrum  data  and  upper 
and  lower  error  bounds  for  each  of  the  runs .   One  program 
using  the  equations  presented  in  Section  IV-A  and  B  computed 
the  elastic  peak  radiation  tail  for  each  run  and  stored  this 
on  cassette  data  files.   Another  program  subtracted  the  radi- 
ation tail  from  the  spectrum  data  and  error  bound  files  for 
each  run  and  created  additional  files  to  store  this  data. 
Appropriate  data  files  were  accessed  by  a  plotting  routine 
and  graphs  were  drawn  using  a  Hewlett-Packard  Model  9872A 
Plotter. 
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The  Hewlett-Packard  program  "REGS"  [Ref.  21]  was  used  as 
a  fitting  routine  to  determine  the  equations  for  the  line 
fit  to  the  inelastic  continuum.   This  program  provides  the 
option  of  specifying  a  linear,  logarithmic,  expotential  or 
polynomial  of  desired  degree. 
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